Unexpected universality in static and dynamic avalanches 
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We find that some equilibrium systems and their non-equilibrium counterparts actually show the 
same jerky response or avalanche behavior on many scales in response to slowly changing external 
conditions. In other words, their static and dynamic avalanches behave statistically the same. This 
suggests that their critical properties are much more generally applicable than previously assumed. 
In this case, systems far from equilibrium may be used to predict equilibrium critical behavior, and 



vice versa. 



PACS numbers: 02.60.Pn, 75.10.Nr, 75.60.Ej 

Avalanche behavior in diverse dynamical systems 
has been extensively studied in the past decade [iHl]. 
In those systems, there are often a large number of 
metastable states. When pushed by an external driving 
field, those systems shift from one metastable state to an- 
other, responding with collective behavior in the form of 
avalanches. A dynamic avalanche is just the rearrange- 
ment of the system configuration, which connects two dif- 
ferent metastable states at two slightly different external 
fields. In experiments, avalanches are often associated 
with crackling noise as measured in acoustic emission and 
Barkhausen noise experiments 2, 3]. So far, avalanche 
behavior in equilibrium systems, i.e. static "avalanche" , 
has rarely been studied due to computational complex- 
ity. With a static "avalanche" we refer to a configura- 
tion rearrangement connecting two different neighboring 
ground states at two slightly different external fields. 

Generally, equilibrium systems are believed to be com- 
pletely different from non-equilibrium ones simply be- 
cause the underlying physics is so different. A natural 
question arises: Do static and dynamic avalanches have 
the same critical behavior? Answering this basic ques- 
tion would be crucial to understand whether there are 
any possible deep connections between equilibrium sys- 
tems and their non-equilibrium counterparts. 

In this letter, we show compelling evidence that static 
and dynamic avalanches have the same critical behavior 
in the zero-temperature random-field Ising model (zt- 
RFIM). This particular model is chosen for two rea- 
sons. First, there is a related highly controversial ques- 
tion in this model, i.e. whether the equilibrium and 
non-equilibrium disorder-induced phase transitions be- 
long to the same universality class. Second, both static 
and dynamic avalanches can be clearly identified and cal- 
culated within this model. We find that all tested uni- 
versal scaling functions and corresponding critical expo- 
nents coincide for static and dynamic avalanches. Our 
findings indicate that generally equilibrium systems and 
their non-equilibrium counterparts may have deep con- 
nections [E, Jj. 

As a prototypical model for disordered magnets, the 
RFIM has been intensely studied [7|. Its Hamiltonian is 



given by 'H = -J '}2<i,j> ^i^j ~ Z^i + where the 
Ising spins Si = ±1 sit on a d-dimensional hypercubic lat- 
tice with periodic boundary conditions. The spins inter- 
act ferromagnetically with their nearest neighbors with 
strength J and experience a uniform external field H and 
a quenched local random field hi. Usually, the local fields 
are chosen from a Gaussian distribution p{h) with mean 
and standard deviation R. R is called the disorder pa- 
rameter. In equilibrium, it is generally believed that in 
d> 2 the transition between the ordered (ferromagnetic) 
and disordered (paramagnetic) phases is continuous and 
controlled by a stable zero-temperature fixed point fsj. 
Therefore, one can set temperature T — Q and tune dis- 
order R to study the equilibrium disorder-induced phase 
transition (DIPT) undergone by the ground-state proper- 
ties. In non-equilibrium, the zt-RFIM has been very suc- 
cessful in explaining dynamic avalanches and crackling 
noise observed in magnets [O-llj. The key result is that 
there is a non-equilibrium DIPT, associated with the hys- 
teretic behavior. Based on the similarities of some critical 
exponents, it has been conjectured that the equilibrium 
and non-equilibrium DIPT may belong to the same uni- 
versality class. But this has been highly controversial due 
to the lack of compelling evidence [a, H, [H, [l^ . Compar- 
ing the critical behavior of static and dynamic avalanches 
is important and necessary to answer this question. 

An avalanche in the RFIM refers to the flip of neigh- 
boring spins during the magnetization process, corre- 
sponding to a jump in the magnetization curve M{H). 
To identify avalanches, we increase the external field 
H from — oo to oo adiabatically slowly, i.e. H is 
kept constant during the propagation of an avalanche. 
Then, a static (dynamic) avalanche connects two nearest 
ground (metastable) states along the equilibrium (non- 
equilibrium) M{H) curve at zero-temperature. In equi- 
librium, the ground state problem of the RFIM can 
be mapped onto the min-cut/max-flow problem of a 
network and solved via the so-called push-relabel algo- 
rithm IJ, [15| . An efficient linear interpolation scheme is 
then used to find steps by narrowing down the H range 
where static avalanches occur [l6i |. In non-equilibrium, 
we use the single-spin-fiip dynamics to calculate the 
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FIG. 1; (Color online) Disorder dependent avalanche behav- 
ior in RFIM. Main panel: Magnetization curves in equilibrium 
(green) and non-equilibrium (red), below, near and above the 
critical disorder Rc. Insets: Cross-sections of 3D systems 
showing all the avalanches (denoted by different colors) oc- 
curring during those magnetization processes. The calcula- 
tion is done on 3D Gaussian zt-RFIM with system size 64^. 
Non-equilibrium: (a) R = 2.0, (b) R = 2.224, (c) R = 2.6. 
Equihbrium: (d) R = 2.25, (e) R = 2.45, (f) R = 2.9. Note 
that RT'^ = 2.16±0.03 and i?"" = 2.28±0.01 for 3D Gaussian 
RFIM 0,[i3|. 



metastable state: each spin flips deterministically when 
its effective local field hf^ = JJ2i^j + + ^ changes 



sign 
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Due to the nearest neighbor interaction, a 
flipped spin will push a neighbor to flip, which in turn 
might push another neighbor, and so on, thereby gener- 
ating an dynamic avalanche. 

To study whether the shape of the random field distri- 
bution would affect the avalanche behavior, we consider 
four different types of p{hys: (1) Gaussian: paih) — 
1 expi-^y, (2) Lorentzian: = j^hj+^k/W'' 

for h £ [-R, R] and else; 



(3) parabolic: pp(/i) 

(4) uniform: p\i{h) — ioi h £ [—R,R] and else. In 
all cases, p{h) is symmetric around h = and the gener- 
alized "width" R will be called the disorder parameter. 

Fig.l shows the M{H) curves and corresponding 
avalanches occurring during the magnetization processes 
at different disorders in both equilibrium (d,e,f ) and non- 
equilibrium (a,b,c). For R < Rc (a,d), most spins tend 
to flip collectively in a system spanning avalanche seen 
as a macroscopic jump in the magnetization curve. For 
R > Rc (c,f), spins tend to flip individually and result 
in many microscopic avalanches and a macroscopically 
smooth magnetization curve. For R ^ Rc (b,e), jumps 
(avalanches) of all sizes are seen in the magnetization 
curve. Qualitatively, we flnd that static and dynamic 
avalanches show similar disorder dependent behavior. 

To quantitatively study the similarity of static and dy- 



namic avalanches, we first study the avalanche size dis- 
tribution integrated over the external field ll| . Near the 
critical disorder Rc, its scaling form can be written as 



(1) 



where S is the avalanche size, i.e. the number of spins 
participating in an avalanche, ± refers to the sign of the 
reduced disorder r = {Rc — R)/R, cr gives the scahng 
of the largest avalanche size S'max kl"^^*^! /? and S 
give the singularities of M{H) near the critical point 
{Hc,Rc)- Here, the critical field He is defined to be 
the field where the slope of M{H) goes to cxd. In non- 
equilibrium, Dint{S, R) for Gaussian p{h) has been stud- 
ied extensively. The critical exponents (r + a/3S) = 
2.03±0.03, a — 0.24±0.02 and the universal scaling func- 
tion &^\X) = e-^-'^^^^'^' (0.021 -t- 0.002X -I- 0.531^2 - 
0.266X'^ -I- 0.261X'*) were obtained from scaling collapses 
of Dint{S,R) at different disorders [11|. 

Fig. 2 (a) shows that for Gaussian, Lorentzian and 
parabolic p(/i)'s, and for both static and dynamic 
avalanches at different disorders, with the same pair of 
critical exponents: (r -|- (7f3S) = 2.03 and a = 0.24, 24 
Dint{S,R) curves collapse onto a single one. The uni- 
versality that the three different p(/i)'s show the same 
avalanche behavior is not a surprise at all. A renormal- 
ization group (RG) analysis has shown that, at least in 
non-equilibrium, what matters is just p"(0), i.e. the sec- 
ond derivative of p{h) at ft, = [lO|. It is easy to check 
that Gaussian, Lorentzian and parabolic p(/i)'s all have 
p"(0) ^ — i?'^. Therefore their universal behaviors agree, 
as expected. Fig. 2(b) shows that for uniform p{h), for 
both static and dynamic avalanches at different disor- 
ders, eight Diat{S,R) curves collapse onto a single one, 
with critical exponents: (t + a(36) = 2.08 ± 0.02 and 
a = 0.52 ± 0.03. Note that: (1) The critical exponents, 
especially a, are significantly different from those of the 
above three kinds of p(/i)'s. (2) The scaling function 
has a significantly different shape from that observed in 
Fig. 2 (a). These two findings are consistent with the RG 
analysis mentioned above because for a uniform p{h), 
p"(0) = is independent of R and is qualitatively dif- 
ferent from the other three distributions. 

The most surprising result about Fig. 2 is that the crit- 
ical exponents and scaling functions for static and dy- 
namic avalanches match for any p{h). This strongly in- 
dicates that the equilibrium and non-equilibrium RFIM 
behave the same near their corresponding critical points. 

To check whether this is just a coincidence, we make 
another independent test by studying the avalanche cor- 
relation function, which measures the probability that a 
distance x between any two fiipping spins occurs in the 
same avalanche ll|. Near the critical disorder Rc, the 



scaling form of the avalanche correlation function inte- 
grated over H can be written as 
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FIG. 2: (Color online) Scaling functions of integrated 
avalanche size distribution. The original Dint{S,R) curves 
for static and dynamic avalanches at different disorders and 
different p{h)'s are calculated in 3D with system size 64"^ 
and are averaged up to 100 initial random-field configura- 
tions. In the legend, the subscripts stand for equilibrium 
(e) or non-equilibrium (n) and the type of p{h): Gaussian 
(G), Lorentzian (L), parabolic (P), and uniform (U). (a) For 
Gaussian, Lorentzian and parabolic p(/i)'s, using critical ex- 
ponents (r-f (7/35) — 2.03 and a = 0.24, 24 curves collapse onto 
each other, up to non-universal critical disorders {R^'^ = 2.16, 
Rf^ = 2.29; i??^ = 1.92, Rf = 2.08; R^ = 4.84, Rf = 5.0) 
and overall scale factors. The thick black curve through 
the collapse is the universal scaling function T>^]}'^{X) of non- 
equilibrium Gaussian RFIM [lj|. (b) For uniform p{h), eight 
curves collapse onto each other, up to non-universal critical 
disorders (ii"^ = 4.64, Rf^ = 4.46) and overall scale factors. 
The collapse yields critical disorders (r -I- cr/35) = 2.08 ± 0.02 
and a — 0.52 ±0.03. Both the critical exponents and the scal- 
ing function are different from those of the Gaussian p{h). 



with V the correlation length exponent and d the dimen- 
sion [l3|. In non-equilibrium, the quantity Gmtix^R) for 
the Gaussian p{h) has been studied extensively, where 
d + p/v = 3.07 ± 0.30 and v = 1.37 ± 0.18 were ob- 
tained from scaling collapses of Gint(a;, R) at different dis- 
orders [111. Here, in Fig. 3(a), we show that for Gaussian, 
Lorentzian and parabolic p(/i)'s, and for both static and 
dynamic avalanches at different disorders, with the same 
pair of critical exponents: d + /3/iy = 3.07 and u = 1.37, 
24 Gint{x,R) curves collapse onto a single one. Fig. 3(b) 
shows that for uniform p{h) and for both static and dy- 
namic avalanches at different disorders, eight Gint{x,R) 
curves collapse onto a single one, with critical exponents: 
d + p/v = 3.0 ± 0.3 and v = 0.8 ± 0.2. On one hand, 
both the critical exponents and the scaling function for 
the uniform p(h) are different from those of the Gaussian 
p{h). On the other hand, both the critical exponents 
and scaling functions for static and dynamic avalanches 
match for any p{h). These results are completely consis- 
tent with what we found in avalanche size distributions. 

Moreover, in a separate work, for the Gaussian p{h), we 
have shown that static and dynamic avalanches have sur- 
prisingly similar spatial structures with the same fractal 
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FIG. 3: (Color online) Scaling functions of integrated 
avalanche correlation function. The original Gint{^, R) curves 
for static and dynamic avalanches at different disorders and 
different p{h)'s are calculated in 3D with system size 64'^ and 
are averaged up to 100 initial random-field configurations, (a) 
For Gaussian, Lorentzian and parabolic p(/i)'s, 24 curves col- 
lapse onto each other, using d + jS/u = 3.07 and v = 1.37 [TTI ]. 
(b) For uniform p{h), eight curves collapse onto each other 
with d + P/u = 3.0± 0.3 and u = 0.8 ± 0.2. Note that those 
collapses are up to the same critical disorders as used in Fig. 2. 



dimensions, anisotropy measures, and associated univer- 
sal scaling functions |18| . More interestingly, we notice 
that the equilibrium and non-equilibrium DIPTs them- 
selves show surprising similarity: (1) They share the 
same no-passing rule: at T — 0, flipped spins can never 
flip back as the magnetic field H is swept monotoni- 
cally I, El]. (2) In mean field theory^ they have the 
same thermodynamic critical exponents [9|, |20| , the same 



avalanche critical exponents [18[ , and the same exponent 
relations [lo|. (3) RG calculations show that the 6 — e ex- 
pansion for the non-equilibrium critical exponents maps 
to all orders in e onto the controversial equilibrium ones: 
The temperature dependence is irrelevant in the equi- 
librium RFIM and the time dependence is irrelevant in 
the zero-temperature non-equilibrium RFIM, leaving us 
with the same starting point for the calculation in both 
cases [l^. All these evidences in favor of universality 
corroborate our findings here. 

To discuss the effect of dynamics on the critical behav- 
ior of avalanches, a general /c-spin-fiip dynamics (with 
fc = 1, 2, • • • , oo) has been introduced 0. It is defined 
such that all the states connected by avalanches are k- 
spin-flip metastable states whose energy cannot be low- 
ered by the flip of any subset of 1, 2, • • • ,k spins (gI. [2l|. 
The case k = 1 just corresponds to the single-spin-flip 
dynamics used in our non-equilibrium calculations. The 
case k = oo corresponds to the ground state evolution dy- 
namics in our equilibrium calculations. It has been found 
that the change of dynamics from k — 1 to k = 2 will not 



alter the critical behavior of the dynamic avalanches [22 1 . 
Together with our flnding, i.e. k = 1 and k = oo give 
the same avalanche behavior, we suggest that avalanches 
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FIG. 4: (Color online) Schematic phase diagram of zt-RFIM 
with fc-spin-fiip dynamics and demagnetization dynamics. 
Dashed Unes stand for the first-order phase transitions occur- 
ring at the critical field H^- Note that for both demagnetiza- 
tion dynamics and ground state evolution (fc = cx3), = 
due to symmetry. (Inset) Schematic M(H) curves associated 
with difi'erent dynamics. 



associated with the whole series of fc-spin-flip dynamics 
(with A; = 1, 2, • • • oo) would have the same critical behav- 
ior. The fc-spin-flip dynamics is quite general but it def- 
initely cannot encompass all kinds of dynamics, e.g. the 
demagnetization dynamics associated with the demagne- 
tization curve, which is obtained by applying an oscillat- 
ing external field with very slowly decreasing amplitude. 
Previous studies for Gaussian p{h) showed that there are 
two very interesting results. First, the avalanches asso- 
ciated with the demagnetization curve are found (within 
numerical error bars) to display the same critical expo- 
nents and scaling functions as the avalanches associated 
with the saturation hysteresis loop (with k = I) [23| . 
Second, the demagnetized state and ground state show 
similarity near their corresponding critical disorders: the 
critical exponents and scaling function associated with 
the M{R) curve coincide [is} . 

Considering all the findings, we suggest that all the 
different dynamics yield the same scaling behavior of 
avalanches — an unexpected universality, see Fig. 4. It 
would be very interesting to numerically test this uni- 
versality in other disordered systems, especially for those 
systems with frustrations where the no-passing rule is 
broken. We suspect that a necessary condition for equi- 
librium and non-equilibrium critical behavior to scale in 
the same way is that the scaling behavior is dominated 
in both cases by a zero-temperature fixed point 



2J. For 



example, for the random-bond Ising model, which has a 
non-trivial finite-temperature fixed point [25|, the equi- 
librium and non-equilibrium critical behavior are differ- 
ent 0. 

We thank James P. Sethna, A. Alan Middleton, Yoshit- 
sugu Oono, Daniel S. Fisher, Gil Refael, Andrei A. Fe- 



[lo: 

[11 
[12: 

[is: 

[14 

[15: 
[16: 

[IT 

[is: 

[19 

[20 
[21 

[22: 

[23: 

[24' 
[25| 



J. P. Sethna, K. A. Dahmen, and C. R. Myers, Nature 
410, 242 (2001). 

G. Durin and S. Zapperi, in The Science of Hysteresis, 
edited by G. Bertotti and I. Mayergoyz (Academic Press, 
2005). 

E. Vives, J. Ortfn, L. Manosa, I. Rafols, R. Perez- 
Magrane, and A. Planes, Phys. Rev. Lett. 72, 1694 
(1994). 

S. Field, J. Witt, F. Nori, and X. Ling, Phys. Rev. Lett. 
74, 1206 (1995). 

A. Maritan, M. Cieplak, M. R. Swift, and J. R. Banavar, 
Phys. Rev. Lett. 72, 946 (1994). 

F. J. Perez- Reche and E. Vives, Phys. Rev. B 70, 214422 
(2004). 

D. P. Belanger and T. Nattermann, in Spin Glasses and 
Random Fields, edited by A. P. Young (World Scientific, 
Singapore, 1998). 

D. S. Fisher, Phys. Rev. Lett. 56, 416 (1986). 

J. P. Sethna, K. A. Dahmen, S. Kartha, J. A. Krumhansl, 

B. W. Roberts, and J. D. Shore, Phys. Rev. Lett. 70, 
3347 (1993). 

K. A. Dahmen and J. P. Sethna, Phys. Rev. B 53, 14872 
(1996). 

O. Perkovic, K. A. Dahmen, and J. P. Sethna, Phys. Rev. 
B 59, 6106 (1999). 

J. P. Sethna, K. A. Dahmen, S. Kartha, J. A. Krumhansl, 

0. Perkovic, B. W. Roberts, and J. D. Shore, Phys. Rev. 
Lett. 72, 947 (1994). 

F. Colaiori, M. J. Alava, G. Durin, A. Magni, and S. Zap- 
peri, Phys. Rev. Lett. 92, 257203 (2004). 
A. A. Middleton and D. S. Fisher, Phys. Rev. B 65, 
134411 (2002). 

1. Dukovski and J. Machta, Phys. Rev. B 67, 014413 
(2003). 

C. Frontera, J. Goicoechea, J. Ortm, and E. Vives, J. 
Comp. Phys. 160, 117 (2000). 

A. K. Hartmann and A. P. Young, Phys. Rev. B 64, 
214419 (2001). 

Y. Liu and K. A. Dahmen, arXiv:cond-mat/0609609v2 

(2006) . 

Y. Liu and K. A. Dahmen, Phys. Rev. E 76, 031106 

(2007) . 

T. Schneider and E. Pytte, Phys. Rev. B 15, 1519 (1977). 

C. M. Newman and D. L. Stein, Phys. Rev. E 60, 5244 
(1999). 

E. Vives, M. L. Rosinberg, and G. Tarjus, Phys. Rev. B 
71, 134424 (2005). 

J. H. Carpenter and K. A. Dahmen, Phys. Rev. B 67, 
020412 (2003). 

D. S. Fisher, private communication. 

K. Hukushima, J. Phys. Soc. Jpn. 69, 631 (2000). 



